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Vacuum energy is a simple model for dark energy driving an accelerated expansion of the universe. 
If the vacuum energy is inhomogeneous in spacetime then it must be interacting. We present the 
general equations for a spacetime-dependent vacuum energy in cosmology, including inhomogeneous 
perturbations. We show how any dark energy cosmology can be described by an interacting vac- 
uum+matter. Different models for the interaction can lead to different behaviour (e.g., sound speed 
for dark energy perturbations) and hence could be distinguished by cosmological observations. As 
an example we present the cosmic microwave microwave background anisotropies and the matter 
power spectrum for two different versions of a generalised Chaplygin gas cosmology. 



I. INTRODUCTION 



Vacuum energy provides a very simple model of dark energy [lMl2j. It is the energy density that remains 
in the absence of any particles, and therefore it remains undiluted by the cosmological expansion. A positive 
vacuum energy can drive an accelerated expansion once the density of ordinary matter or radiation becomes 
sub-dominant. Unlike other models of dark energy it does not necessarily introduce any new dynamical degrees 
of freedom. 

We define a vacuum energy, V, to have an energy-momentum tensor proportional to the metric 

tji = -Vtf . (1) 

By comparison with the energy-momentum tensor of a perfect fluid 

T£ = Pg? + ( P + P)u»u v , (2) 

we identify the vacuum energy density and pressure with p = —P = V, but since there is no particle flow then 
the four- velocity of the vacuum, is undefined. 

A vacuum energy that is homogeneous throughout spacetime, V^V = 0, is equivalent to a cosmological 
constant in Einstein gravity, A = 8ttGnV. The discrepancy between the value of the energy density required by 
current observations and the typical energy scales predicted by particle physics is the long-standing cosmological 
constant problem 

We will consider the possibility of a time and/or space dependent vacuum energy. From Eq. ([T]) we have 

V M T£ = Q u . (3) 

where the energy flow is given by 

Q y = -V v V. (4) 

We can therefore identify an inhomogeneous vacuum, V^V ^ 0, with an interacting vacuum, Q v ^ |l[. The 
conservation of the total energy-momentum (including matter fields and the vacuum energy) in general relativity 

V, (T£ + f?) = , (5) 

implies that the vacuum transfers energy-momentum to or from the matter fields 

V M T^ = -Q v . (6) 

Note that the energy density and four- velocity of a fluid can be identified with the eigenvalue and eigenvector 
of the energy- momentum tensor ^ : 

Tjlu" = -pu^ . (7) 

Because the vacuum energy-momentum tensor ([T|) is proportional to the metric tensor, any four- velocity, u M , is 
an eigenvector 

fy u = -Vu" V , (8) 
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and all observers see the same vacuum energy density, V, i.e., the vacuum energy is boost invariant. 

Although the vacuum does not have a unique four- velocity, we can use the energy flow, Q„, to define a 
preferred unit four- vector in an inhomogeneous vacuum [l| 

normalised such that u^ii^ = ±1 for a spacelike or timelike flow. Note however that the defines a potential 
flow, i.e, with vanishing vorticity. 

In this paper we will consider vacuum energy which may be inhomogeneous in spacctime, interacting with fields 
or fluids without necessarily invoking additional degrees of freedom. We show that any spatially homogeneous 
dark energy cosmology can be decomposed into an interacting vacuum+matter cosmology. By lifting this 
spatially-homogeneous solution to a covariant interaction one can study inhomogeneous perturbations which 
obey coupled first-order equations of motions for the matter density and velocity. As an example we consider 
perturbations of a generalised Chaplygin gas cosmology, decomposed into an interacting vacuum+matter. 



II. VACUUM COSMOLOGY 



A. FRW background 

The symmetries of a spatially homogeneous and isotropic Friedmann-Robertson- Walker (FRW) metric, with 
scale factor a(t) and Hubble rate H = a/a, require the vacuum to be spatially homogeneous and isotropic too, 
hence V = V(t). In this case the vacuum and matter are both homogeneous on spatial hypersurfaces orthogonal 
to the matter four-velocity, = (1, 0, 0, 0). Note that the energy flow, u^, and matter velocity, u M , necessarily 
coincide in FRW cosmology due to the assumption of isotropy. 

The Friedmann constraint equation requires 

& = ^(p + V)-^, (10) 

where K determines the spatial curvature. The continuity equations for matter fields and vacuum are 

p + 3H(p + P) = -Q, (11) 
V = Q. (12) 

The vacuum energy, V, is undiluted by the cosmological expansion, but can have a time-dependent density in 
the presence of a non-zero energy transfer, Q ^ 0, where 

Q = -u v Q v . (13) 

There have been many attempts to describe the present acceleration as due to a time-dependent vacuum 
energy [||-[l2j]. However, there is little to be learnt from simply assuming an arbitrary time-dependent vacuum 
to obtain the desired cosmological solution. Ideally one should have a physical model from which one can derive 
a time-dependent solution and study other physical effects. For example, vacuum fluctuations of free fields can 
support an averaged density proportional to the fourth-power of the Hubble expansion, V ex H A [14j . Such a 
vacuum e nerg y would not in itself support an accelerated expansion, but other forms such as V oc H have been 
proposed [15( better able to match the observational data [161420] . 

It is not obvious how such time-dependent vacuum models can be compared with observations in an inhomo- 
geneous universe. We will argue that it is possible to give a consistent description of vacuum dynamics, and in 
particular the relativistic equations of motion for inhomogeneous perturbations, given a covariant, physical pre- 
scription for the local vacuum energy or, equivalently, the vacuum energy transfer 4-vector, = — V M V. One 
should then be able to subject vacuum models to observational constraints, even in the absence of a Lagrangian 
derivation or microphysical description. 



B. Linear perturbations 



Let us consider inhomogeneous linear, scalar perturbations where the energy and pressure of matter is given 
by p(t) + fipfy, x l ) and P(t) + SP(t, x l ), and the four-velocity of matter is given by 

u»=[l-<l>,a- 1 d i v] , Ufl = [-l-<t>,di6] . (14) 

where we define d l v = a(dx l /dt) and 6 = a(v + B). Once we allow for deviations from homogeneity in the 
matter and metric, we should also allow for inhomogeneity in an interacting vacuum, V(t) + 5V(t,x l ). As 
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remarked earlier, the vacuum has an energy density and pressure, but no unique velocity. In particular the 
momentum of the vacuum vanishes in any frame, (p + P)9 = 0. On the other hand the energy flow u defined 
in Eq. ©, can be written in analogy with the fluid velocity (fT4")l as 

fi" = [1 - , a-'d'v] , fi M = [-1-0, di§\ . (15) 

where from Eq. (J9j> we identify 9 = —SV/V. 

Perturbations about a spatially flat (K — 0) FRW metric plri^ij are described by the line element 



ds 2 = -(1 + 2c/))dt 2 + 2ad l Bdtdx i + a 2 [(1 - 2ip)Sij + 2d l d 3 E] dx l dx 2 . (16) 
Following [2l], [12, [HJ , we decompose the energy-flow along and orthogonal to the fluid velocity, 

Q» = + U , (17) 

where = 0, so that we have 

Q ll = [-Q{l + cl>)-8Q,d i (f + Q6)]. (18) 
The energy continuity equations for matter and vacuum become 

Sp + 3H{Sp + SP)-3{p + P)iJj + {p + P) — (0 + a 2 E-aB) = -SQ - Q<j> , (19) 

SV = 5Q + Q<j). (20) 

while the momentum conservation becomes 

(p + P)8-3c 2 s H(p + P)9 + {p + P)4> + 6P = -f + c 2 s Q8, 

SV = f + QO. (21) 

where the adiabatic sound speed c 2 = P '/ p. Note that the vacuum momentum conservation equation becomes 
a constraint equation which requires that the vacuum pressure gradient is balanced by the force 

Vi(-V) = Vi(/ + Q0). (22) 

This determines the equal and opposite force exerted by the vacuum on the matter: 

- f = SV + V6 . (23) 

i.e., the fluid element feels the gradient of the vacuum potential energy. 

Note that the perturbations of a fluid coupled to the vacuum with P = —p has no additional degrees of 
freedom, in contrast to a dark energy fluid with Px =/= —px- Using the vacuum energy and momentum 
conservation equations to eliminate SQ and / we obtain 

Sp + 3H(Sp + SP) - 3(p + P)ip + (p + P)\ (9 + a 2 E - s) - SV , (24) 
(p + P)0 - 3c 2 s H{p + P)6 + {p + P)cp + SP = SV + (1 + c 2 s )V6 . (25) 



C. Gauge invariant perturbations 

It is well known that metric and matter perturbations can be gauge-dependent under a first-order gauge 
transformation, such as t — > t + St(t,x l ). The fluid density and pressure transform as Sp — >• Sp — pSt and 
SP — > SP — PSt [24]]. Similarly the vacuum perturbation transforms as SV — > SV — QSt and SQ — > SQ — QSt. 
The fluid 3-momentum transforms as 6 — > 9 + St and the energy flow transforms similarly as 9 — > 9 + St. 

We can construct gauge-invariant perturbations by specifying quantities on a particular physical reference 
frame (24[. For example, the vacuum perturbation on hypersurfaces orthogonal to the energy transfer, Q M , can 
be shown to vanish identically: 

AV com = SV + V6 = , (26) 

since from Eq. (JSJ and (fT5")) we have 9 = —SV/V. This simply reflects that the energy flow is the gradient of 
the vacuum energy and therefore the orthogonal hypersurfaces are uniform vacuum energy hypersurfaces by 
construction. 



(b) 




FIG. 1: (a) In an FRW cosmology the homogeneous spatial hypersurfaces are orthogonal to both the fluid 4-velocity, 
u M , and the vacuum energy flow, Q M . (b) In an inhomogeneous cosmology the spatial hypersurfaces orthogonal to the 
fluid 4-velocity (light orange) and the vacuum energy flow (dark red) do not necessarily coincide. 



On the other hand the vacuum perturbation on hypersurfaces orthogonal to the matter 4-velocity, w M , (the 
comoving vacuum perturbation) is given by 

SVcam = 5V + V8 = -f . (27) 

This is in general non-zero, i.e., the vacuum may be spatially inhomogeneous in the comoving-orthogonal gauge. 
For example, the Poisson equation for the Newtonian metric potential is given by 

V 2 $ = AttG {6p coia + 6V com ) , (28) 

i.e., the Newtonian metric potential is sourced by both the matter and vacuum perturbations, where 

Spcom = Sp + p6. (29) 

Note that we can write the comoving vacuum density perturbation (|27[) as 

5V com = V(9-9) . (30) 

Therefore, if the energy flow follows the fluid four- velocity, u M = u M , then we have 6 — 6 and the vacuum is 
spatially homogeneous on comoving-orthogonal hypersurfaces, SV com = 0. 

Another gauge invariant expression for the vacuum density perturbation is the dimcnsionless vacuum pertur- 
bation on uniform-fluid density hypersurfaces, which describes a relative density perturbation 
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If, for example, the vacuum energy is a function of the local matter density, V — V(p), then the relative density 
perturbation must vanish and the vacuum is spatially homogeneous on uniform-density hypersurfaces, S — 0. 

The total non-adiabatic pressure perturbation due to any intrinsic non-adiabatic pressure of the matter and 
the relative entropy perturbation between the vacuum and matter is then 

<5P„ad = SP 

= SP 

where the adiabatic sound speed for matter is c 2 s = P/p. This vanishes for adiabatic matter perturbations, 
5P = c 2 Sp, and adiabatic vacuum fluctuations, S = 0, or a non-interacting vacuum, Q = 0. 



-SV- 



-c 2 Jp- 



P-V 



(Sp + SV) , 



{l + c 2 s )Q[Q + 3H{p + P)} 
9H 2 (p + P) 



S, 



(32) 



III. DECOMPOSED GENERALISED CHAPLYGIN GAS 



As an example of how an inhomogeneous vacuum energy might be used to describe the present accelerated 
expansion of our Universe we will show how one widely-studied dark energy model, the Chaplygin gas, can 
be re-interpreted as an interacting vacuum+matter cosmology, and how this re-interpretation can motivate 
different possible behaviour for density perturbations. 

Any dark energy fluid energy-momentum tensor |(2J) with density pdc can be described by pressureless matter, 
with density p m and velocity uf n = m m , interacting with the vacuum, V, such that pdc = Pm + V The 
corresponding matter and vacuum densities are given by 

Pm, = Pdc + Pdc , V = -Pdc • (33) 

while the energy flow is — V^Pdc- In an FRW cosmology this corresponds to Q = —Pdc- One might choose 
to decompose a dark energy model Pde(a) into any two interacting barotropic fluids such that pdc = Pi + Pi-, 
but this would double the degrees of freedom in the model unless one of these two "fluids" is the vacuum. 
The generalised Chaplygin gas, defined by the barotropic equation of state [2^, [26[ 

P s c g = -Ap g « g . (34) 

This leads to a solution for the density in an FRW cosmology 

/ , n\ 1/(1+") 

PsCs =(A + Ba-^) \ (35) 

This has the simple limiting behaviour p g c g oc a~ 3 as a — > and p g c g — > A 1 /( 1+a ' > as a — > +oo, therefore this has 
been proposed as a unified dark matter model. However such models are strongly constrained by observations 
since the barotropic equation of state defines a sound speed for matter perturbations which only reproduces the 
successful ACDM model when a -> [13, HI- 

The decomposition (|33[) into pressureless matter interacting with the vacuum has previously been considered 
for the generalised Chaplygin gas by Bento et al [29]. In this case we have the FRW solution 

(\ —a/ 
A + Ba-'^ 1+a n , (36) 

and hence 



A=( Pm + V) a V. (37) 
The form of the FRW solution suggests a simple interaction 

In the matter or vacuum dominated limits this reduces to an interaction of the form Q oc Hp m or Q oc HV 
studied, for example, by Barrow and Clifton [30(. 

It is intriguing to note that the FRW solution can be defined in terms of an interaction (|38p with a single 
dimensionless parameter, a, whereas when defined in terms of an equation of state (|34j) its definition requires 
both a and the dimensional constant A which determines the late-time cosmological constant. In the interaction 
model, A (and therefore the late-time cosmological constant) emerges as an integration constant dependent on 
initial conditions. 
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However, to study inhomogeneous perturbations in the decomposed model we must "lift" the explicitly time- 
dependent FRW solution to a covariant model for the interaction. We have at least two choices. In either case 
one can calculate the speed of sound in the combined interacting vacuum+matter using the usual definition 
[3TM33L l39l ] and the decomposed density and pressure (l33l) 



2 _ fSP d0 \ ( -sv 



-u^j ■ (39) 



A. Barotropic model 

Firstly one could require that the local vacuum energy is a function of the local matter density, V = V(p m ). 
If this applies to inhomogeneous perturbations of the interacting vacuum+matter as well as the background 
then we require 

SV = —5 Pm . (40) 

Pm 



This implies that the vacuum perturbations are adiabatic and S — in Eq. pip . On the other hand the 
comoving vacuum energy, (|27[). is non-zero 

V 

SV com = —8pm.com ■ (41) 
Pm 

Thus one needs to consistently include the vacuum inhomogeneities as well as matter inhomogeneities, for 
example in the Poission equation (|28p. when 7^0. 

Using the adiabatic condition (l40l) we obtain the dark energy sound speed ((39)) 

cf nt = 5^ ■ (42) 

PgCg 

We see therefore the the sound speed for the interacting vacuum+fiuid is the same as the adiabatic sound speed 
of the original barotropic Chaplygin gas (|34p . Perturbations of the the interacting vacuum+matter respect the 
same barotropic equation of state as the original fluid (|34|) and therefore the observational predictions for the 
cosmic microwave background (CMB) anisotropies, for example, are exactly the same as the original fluid model 

& 

CMB temperature anisotropies as well as the power spectrum for the interacting vacuum+matter are shown 
in figures 2 and 3. We see that only very small values of the dimensionless parameter \a\ < 10~ 4 are allowed 
without introducing unacceptably large oscillations in the matter power spectrum [27], , forcing the model to 
be extremely close to the standard ACDM cosmology. 



B. Geodesic flow 



The interacting vacuum+matter model allows a different behaviour for inhomogeneous perturbations than 
the barotropic fluid. Suppose that the energy flow, Q^, is along the matter 4- velocity, u^. The force exert 
exerted by the vacuum on the matter, — / M in Eq. (1171) . is zero and matter 4- velocity is a geodesic flow. 

In this case the comoving vacuum perturbation, 5V com in Eq. (|27p . is zero, i.e., the vacuum is spatially 
homogeneous on hypersurfaces orthogonal to the matter 4-velocity. This means that the comoving pressure 
perturbation for the interacting vacuum+matter is zero and hence the sound speed (|39p is zero, cf nt = [35l.l36j|. 

The evolution of density perturbations is therefore much closer to that in a standard ACDM cosmology with 
zero sound speed than the original Chaplygin gas model with a / 0. CMB temperature anisotropies as well 
as the matter power spectrum are shown in figures 4 and 5. Much larger values of a ~ 0.1 may be compatible 
with the data in this case [34| . 

Note however that this geodesic model for the interaction does allow a gauge-invariant vacuum entropy 
perturbation (|3"Tjl 

S = -, — 5p mtCom . (43) 

Pm 

However due to the Poisson constraint equation (1281) the comoving density perturbation vanishes on large scales 
and hence this is compatible with adiabatic initial conditions for the interacting vacuum+matter at early times. 

The requirement that = Qu^ does impose an important physical restriction on the matter 4-velocity 
since the energy flow must then be irrotational, i.e., it M cx V^V. This would have important consquences 
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FIG. 2: Cosmic microwave background angular power spectrum for the generalised Chaplygin gas with barotropic 
equation of state and an adiabatic sound speed [341 ]. 



10 H 



10" 



10' 



10' 



& 10° 



g 10" 

CL 

10" 



10"% 



10 



10 



a=0.05 

■ a=0.01 
a=0.001 
a=0.0001 
a=0.00001 

■ a=0.000001 




10" 



10" 



10"* 

k[h Mpc -1 ] 



10" 



10 u 



FIG. 3: Power spectrum for baryons plus generalised Chaplygin gas with barotropic equation of state and an adiabatic 
sound speed [34| . 
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for non-linear structure formation, possibly disastrous if gravitationally collapsed dark matter halos could not 
be supported by rotation. Hence we only expect geodesic matter interacting with the vacuum to be a viable 
model for the initial stages of structure formation and we would need a microphysical model for gravitationally 
collapsed halos. 



IV. DISCUSSION 



Many different explanations have been proposed for the present-day acceleration of the Universe [2] . Perhaps 
the most common is a self-interacting scalar held, ip, whose self-interaction potential energy could dominate the 
present energy density. Different models are defined by the chosen functional form for the potential energy, V(<p) 
[37L l38l|. or the kinetic energy as a function of the field gradient, X = (Vy ?) 2 l39l. Eoj . Alternative explanations 
include fluid models, specified by a barotropic equation of state, P(p) [25|, [4j], |42| , some of which could provide 
unified dark matter models capable of explaining acceleration on large scales and galactic dynamics (as dark 
matter) on much smaller scales. Acceleration requires an exotic equation of state with negative pressure, and yet 
the sound speed must remain real in order to avoid instabilities with respect to inhomogeneous perturbations. 
More sophisticated models have been proposed which allow for interactions between the different fields or fluid 
components, e.g., coupled quintessence models p3-48]. Such models are motivated by astronomical observations 
0, providing increasingly detailed constraints on the model parameters, but the models themselves lack a 
persuasive underpinning physical motivation. One might therefore simply consider the simplest possible model, 
or parameterisation, compatible with the data. 

In this paper we have considered vacuum energy as a source of spacetime curvature in Einstein gravity. 
A homogeneous vacuum in Einstein gravity is equivalent to a cosmological constant, whereas an inhomoge- 
neous (time- or space-dependent) vacuum implies an interacting vacuum. Simply specifying a particular time- 
dependence for the vacuum energy in an FRW cosmology is unsatisfactory if the vacuum energy is introduced 
solely to produce a particular time-dependence of the cosmological expansion. Different physical models for 
the origin of the vacuum energy, or its interaction with other matter fields, will lead to different cosmological 
behaviour. Even models which yield the same FRW background solutions [§3, [H| may be distinguished for 
instance by the predictions for CMB anisotropics or the matter or galaxy power spectrum. 

As an example, we have shown how the generalised Chaplygin gas cosmology can be re-derived as a solution 
to an interacting mattcr+vacuum model. The interaction can be defined by a single dimensionless constant and 
the late-time constant vacuum energy (which appears as a dimensional parameter in the original Chaplygin gas 
model) can instead be derived as a constant of integration in the matter+vacuum model. We have shown the 
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FIG. 5: Power spectrum for baryons plus the decomposed generalised Chaplygin gas with geodesic matter and zero 
sound speed [34| . 



CMB and matter power spectrum predictions for two models for interacting vacuum+matter which yield the 
same background solution as the generalised Chaplygin gas, but give very different observational predictions 

Another familiar example, which we have not discussed here, would be a quintessence model with a self- 
interacting scalar field, ip. The self-interaction potential of the field, V(ip), provides a vacuum energy density 
interacting with the kinetic energy of the scalar field. It is well known that a canonical scalar field has a sound 
speed equal to unity and one might thereby hope to distinguish it from a barotropic fluid model where the 
sound speed is necessarily equal to the adiabatic sound speed. By considering a broader class of interacting 
vacuum models one can study models with a range of different possible sound speeds. Our equations enable us 
to consider vacuum energy interacting with other forms of matter, including pressureless matter or radiation. 
We can therefore consider vacuum energy models which do not introduce any degrees of freedom beyond those 
already present in the cosmology, e.g., unified dark energy models. 

In our work we have considered only linear perturbations but non-linear evolution may provide further 
observational constraints on different models. Often it is assumed that the vacuum energy remains unperturbed, 
or negligible, during collapse, but some cases, such as the interacting vacuum+matter with geodesic 4- velocity 
and zero sound speed may be amenable to a study of non-linear collapse [49| and we hope to study this further 
in future. 
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jara, for their hospitality. DW is supported by STFC grant ST/H002774/I. JD-S is supported by CONACYT 
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